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Abstract. We report 51 V zero-field NMR of manganese vanadate spinel of MnV204, 
together with both ac and dc magnetization measurements. The field and temperature 
dependence of ac susceptibilities show a reentrant-spin-glass-like behavior below the 
ferrimagnetic(FEM) ordering temperature. The zero-field NMR spectrum consists of 
multiple lines ranging from 240 MHz to 320 MHz. Its temperature dependence reveals 
that the ground state is given by the simultaneous formation of a long-range FEM 
order and a short-range order component. We attribute the spin-glass-likc anomalies 
to freezing and fluctuations of the short-range ordered state caused by the competition 
between spin and orbital ordering of the V site. 
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1. Introduction 

For several decades, transition-metal spinels of AB 2 X 4 type have been a subject of 
active research because of the intriguing nature of the underlying physics pQ. The 
many unusual phenomena found in the spinel compounds are due primarily to the 
corner-shared tetrahedral network of the B cations which resembles the geometrically 
frustrated pyrochlore lattice. The magnetic properties of spinels rely heavily on whether 
the B cations have orbital degrees of freedom, and whether the A cations are magnetic 
or nonmagnetic. 

Among the vast number of spinels, manganese vanadate MnV204 (2, |3] is one of the 
fascinating but very complex compounds to understand. This is because (i) the V 3+ ion 
(3d 2 , S — 1) possesses orbital degeneracy in its t 2g orbital; the orbital ordering being a 
common feature found in other insulating vanadates (A = Zn jl], Mg [5], Cd [6]); and 
(ii) both the A and B sites are magnetic as in manganese chromite, MnCr 2 04 [£]• 

One may expect that MnV204 shares common features with MnCr204 as well as 
with other vanadium spinels. Indeed, MnV204 undergoes long-range FEM ordering at 
56 K and the spin configuration of the B sites changes from collinear to non-collinear 
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triangular type at low temperature (T), as found in MnCr 2 04 [8]. Like other vanadate 
spinels, MHV2O4 also undergoes two successive phase transitions - orbital order and 
magnetic long-range order. Interestingly, in the case of M11V2O4, the cubic-to-tetragonal 
structural transition related to the orbital order takes place at T$ = 53 K, just below 
the FEM transition temperature Tq = 56 K [3j. Here the question arises about the 
combined effect of the two features (i) and (ii) on a ground state. However, this issue 
has not yet been thoroughly addressed. 

In this paper, we report the observation of the complex ground state in MnV204 
through ac/dc magnetization and zero-field (ZF) 51 V NMR measurements. Our data 
reveal that there exists a short-range ordered state on top of the long range FEM order, 
leading to a spin-glass-like behavior. They are discussed in terms of the competing spin 
and orbital ordering of the V 3+ ion. 

2. Sample preparation and experimental details 

Polycrystalline powder sample of MnV204 was prepared by standard solid-state 
reaction. Stoichiometric mixtures of MnO and V2O3 were ground together and pressed 
into pellet, the pellet was placed into an evacuated quartz tube (~ 1CT 5 torr) and 
fired for 40 hours at 950 °C. Samples were characterized by X-ray powder diffraction. 
Both dc and ac susceptibility measurements were performed as a function of both field 
and temperature using SQUID magnetometer (Quantum Design MPMS). An applied 
frequency was in the range of 10 _1 — 10 3 Hz with a driving field magnitude of 5 
Oe. NMR experiments were carried out at zero external field using a coherent pulsed 
spectrometer capable of computer-controlled tuning of a probe which is calibrated over 
a wide frequency range. In preliminary measurements using an untuned NMR probe 
(Q = 1), we detected two groups of signals at near 280 MHz and 560 MHz at 4.2 K. 
These were assigned to the NMR signals of V 3+ and Mn 2+ ions, respectively. A probe 
with tank circuit was subsequently employed to focus on the stronger signal of the V 
site. The 51 V NMR spectra were obtained by integrating averaged spin echo signals as 
the frequency was swept through the resonance line. 

3. Results and discussion 

3.1. Magnetic susceptibilities 

Figure 1 shows the temperature dependence of the dc susceptibility under an applied 
field of 300 Oe together with the ac susceptibility at a frequency of 1 Hz. The field- 
cooled (FC) and zero-field-cooled (ZFC) curves display a pronounced bifurcation below 
the FEM ordering temperature of Tq- As Tq is approached from above, the FC 
curve increases steeply and then shows a dip around the structural phase transition 
temperature at T$ and finally increases monotonically. In contrast, the ZFC curve 
exhibits a maximum and then a monotonic decrease upon cooling from Tq- Our results 
are in a good agreement with Ref. [3j. Noticeably, one observes a substantial difference 
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Figure 1. (Online color) Temperature dependence of the ac (circle symbols) and dc 
(solid lines) susceptibilities of M11V2O4. For clarity, the real and imaginary part of the 
ac susceptibilities are divided by 2 and multiplied by 6, respectively. 




Figure 2. (Color online) Field dependence of the real (a) and imaginary (b) parts 
of the ac susceptibility of M11V2O4 as a function of temperature. A field is driven at 
H ac = 5 Oc and v = 100 Hz. 



between the ac and ZFC dc susceptibilities below Tq, implying the presence of reentrant- 
spin-glass (RSG)-like anomalies. 

In figure 2 we present the temperature dependence of the ac susceptibilities, 
x{H a , T) for different static field H a applied parallel to the ac driving field. With 
decreasing temperature the real part of x{H a ,T) exhibits a sharp peak around Tc and 
then a round maximum and finally falls off upon further cooling. With increasing H a , 
the latter maximum (designated by Th) decreases in both amplitude and temperature. 
This so-called Hopkinson maximum is not critical in origin, but arises from processes 
associated with the regular /technical contributions to the susceptibility (for example, 
domain wall motion and coherent rotation) [9] . The former sharp peak develops to two 




Figure 3. (a) The critical field amplitude versus the internal field (open square) 
is presented together with the reduced temperature versus the internal field (full 
triangle) on a double-logarithmic scale, (b) A double-logarithmic plot of the critical 
field temperature versus the reduced temperature. The solid and dotted lines are a fit 
toEqs. (l)-(3). 



peaks (denoted by T m and T$, respectively) with increasing dc field. The peak T m is 
suppressed in amplitude and shifts to higher temperature with an increase in the dc field. 
This is due to critical fluctuations accompanying a transition from a paramagnetic to a 
FEM state [10]. In contrast, the peak T$ decreases in amplitude without showing any 
shift in temperature. This might be associated with the structural phase transition [3], 
leading to a pinning of domain walls. The imaginary part of x{H a , T) shows a sharp 
peak around Tc, which exhibits a similar behavior as the peak T$. In addition, the loss 
peak of x"(H a ,T) shows up around 51 K, which is suppressed in both temperature and 
amplitude as T a increases. 

Since the peak T m is governed by critical fluctuations, it is expected to follow the 
static scaling law. The field and temperature dependence of the critical peak is related 
to the standard critical exponents [Tlj [T2j IT5] by 



X^T^cxH] 15 - 1 (1) 
t m oc H^' 1 (2) 

X (Hi,T m )<xtf, (3) 

where Hi is the internal field given by H a = Hi — NM where N is the demagnetizing 
factor and M is the magnetization. In our case, N is estimated by the slope of the low 
field shearing curves around Tc (not shown here). t m is the reduced temperature defined 
as t m = (T m — T c )/T c . We determine the temperature of T c = 56.58 K from the zero 
field data. The critical field amplitude x(Hi,T m ) is directly obtained from Fig. 2(a). In 
Fig. 3 we provide a double-logarithmic plot of x(Hi,T m ) versus Hi, t m versus Hi, and 
x(Hi,T m ) versus t m , respectively. A close look at the data reveals a small curvature at 
low fields. This is typical for the system with exchange coupling strength disorder [TO] . 
A least squares fit of the data between the restricted fields yield the exponent values of 
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Figure 4. (Color online) Frequency dependence of the real (a) and imaginary (b) 
parts of the ac susceptibility as a function of temperature. Inset: Frequency shift of 
the inflection point of the real part of the susceptibility versus temperature. 



5 = 4.47 ± 0.05, 7 + (3 = 0.95 ± 0.04 and 7 = 0.59 ± 0.03. The 5 and (3 values are close 
to the prediction of a 3D Heisenberg model p2] . However, the 7 value strongly deviates 
from it. Although we cannot exclude the uncertainties in evaluating exponents, it might 
be ascribed to a first-order-like transition due to strong spin-orbital coupling as well as 
to a competition between frustration and a structural phase transition (see below). 

In figure 4 the temperature dependence of the x'(T) and ~x!'(T) is presented as a 
function of frequency in the range of 10 _1 — 10 3 Hz. Both x'(T) an d x"{T) exhibit 
pronounced frequency dependence between 20 K and Tq- Overall, with increasing 
frequency both x'(T) and x"(T) shifts to higher temperature. This is unexpected for a 
long-range ordered FEM state, implying that a spin freezing phenomenon takes place 
in the FEM ordered background. 

The shift of the ac susceptibility as a function of the frequency is estimated by 
taking the inflection point of x'(T) [see the inset of figure 4(a)]. From the maximal shift 
ATi we obtain the value of tp = AT/TA( logu) ~ 0.007. For a spin glass the value is 
typically of the order of 10~ 2 , while for a superparamagnet it lies between 10 _1 — 10 -2 
[14] . Using this criterion, we conclude that the studied compound belongs to a class of 
a spin glass. Further, we fit the experimental data in terms of the Vogel-Fulcher law 
uj = u exp[—E a /k B (Ti — T )] to figure out the nature of the spin freezing process. It 
yields the values of ujq w 10 33 Hz, E a pa 287 K, and To ~ 33.5 K. The unrealistic large 
value of u suggests that the spin freezing is not of simply cooperative but governed by 
an intricate process. 



3.2. NMR measurements 

51 V zero- field NMR spectrum and its temperature evolution are shown in figure 5(a) for 
slow cooling (SC), with the rate of ~1 K/min, and in figure 5(b) for rapid cooling (RC), 
with the rate of 60 K/min. The Boltzmann correction for the signal intensities has been 
made by multiplying each spectrum by T. The spectrum has a complex structure and 
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spans a wide range of frequencies from 240 to 320 MHz. The resonance frequency in zero 
field is given as •JnHm where 7^ (= 11.193 MHz/T) is the nuclear gyromagnetic ratio. 
The hyperfme field ifhf is dominated by the core polarization of the inner s-electrons 
by the outer unpaired <i-electrons. For the V 3+ ion (3d 2 , S = 1), the estimated ifhf due 
to pure Fermi contact is about ~ 25 T [15J which falls into ~ 280 MHz range. The 
contributions from transferred, dipolar, or orbital hyperfme fields are negligible. 

One can infer from figure 5 that the 51 V NMR line is composed of several overlapping 
structures, whose peaks are somewhat resolved. This is not surprising given the low 
crystal symmetry [16J. One possibility is the quadrupole interactions split the 51 V 
(/ = 7/2) degenerate levels producing 21 = 7 transitions. However we discard this 
model since the quadrupole coupling required to represent the spectrum would be too 
large for the V nucleus. Hence, we attribute the structures to V sites belonging to 
several magnetically differentiated sites. 

The temperature dependence of the integrated intensity multiplied by T of the 
spectrum is shown in figure 5 (c) for both SC and RC. The signal for the RC case can 
be observed even above Tq while the signal for SC becomes very weak near 40 K. The 
loss of NMR signal intensity at high T is due to the fast spin fluctuations i.e., shortening 
of T2 associated with magnetic ordering. The different temperature dependence of the 
signal intensity for different cooling rate points to the effects of disorder which influences 
the freezing rate of the short-range correlated spins, a conclusion also deduced previously 
from the magnetization results. At low T < 20 K, however, the signal intensity and 
shape are the same regardless of the cooling rate implying that the spin ground state 
of this compound is not affected by the thermal history. This is consistent with the 
ac susceptibility results, which show no frequency dependence at the same temperature 
range. 

By examining the evolution of the spectrum in figure 5 (a) and (b), we take note 
that the NMR peak near 265 MHz is quite distinguishable from the other lines, and 
survives after the disappearance of the other lines at high temperature. Also we see 
that the shift of the NMR line is very small, if any, up to temperature near Tq- This is 
surprising since the NMR resonance frequency usually follows the bulk magnetization 
which often shows gradual second order transition behavior. Therefore, the NMR data 
suggest that the magnetic transition at Tc must be of first order. 

In both SC and RC cases, the stable configuration is established below 15 K as 
shown in figure 5. Here we show that the intriguing ground state is expected in the 
presence of several competing interactions. Spinels contain two different exchange 
coupling constants: Jab and Jbb between the A and B and the B and B sites, 
respectively. For a cubic phase, a magnetic configuration relies on the parameter u 
defined as [T7j 

u — ^^ BB ^ B (4) 
3</abSa 

where S\ and Sb are spin magnitudes at the respective A and B sites. The relative 
strength of two interactions determines a ground state. For example, the ground state 
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Figure 5. (Online color) (a) Evolution of the NMR spectrum with increasing T 
in ZF for slow cooling condition. The Boltzmann correction for intensity was taken 
into account by scaling each spectrum by T. (b) For rapid cooling, the evolution of 
the spectral shape is the same except the observation of the signal at much higher 
temperature, (c) Total integrated intensity of the NMR signal in zero field was plotted 
against temperature for slow and rapid cooling. Note that, for rapid cooling, the signal 
was detected even above Tq- 
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of chromites [Co,Mn]Cr204 with u — 1.5 — 2.0 shows a coexistence of FEM long-range 
order and spiral short-range order that causes the RSG-like behavior [Sj. The M11V2O4 
system seems to lie in the same parameter range as suggested by the fact that the 
RSG-like behavior starts to appear at Tc of the cubic phase (see figure 1). 

In [Co,Mn]Cr 2 4 the appearance of the short-range order is ascribed to the residual 
magnetic frustration in the B site, which is non- vanishing when Jab is smaller than Jbb- 
In our case, however, this scenario seems to be less plausible since the frustration could 
be relived through the orbital ordering of the V 3+ ions. Rather, the short-range character 
in MnV20"4 may be due to the specific orbital ordering that occurs at Tg. In the cubic 
phase (T > Tg) the t 2g orbitals are degenerate and thus the orbital degrees of freedom 
are irrelevant. The V sublattice aligns ferromagnetically due to the antiferromagnetic 
couplings between V 3+ and Mn 2+ ions. In contrast, in the tetragonal phase (T < Tg), a 
structural phase transition takes place to one whose orbital configuration is intimately 
coupled to a spin configuration via a orbital-spin coupling [3 J . Accordingly, an antiferro- 
orbital arrangement of the V 3+ ions is preferred due to the dominant ferromagnetic 
interaction between V 3+ ions. 

This staggered-type orbital configuration of MhV^O^ is contrasted to the orbital 
patterns found in vanadium spinels (A=Zn,Mg,Cd), which are the nonmagnetic analogue 
of MnV 2 4 in the A site [TBI EH 120] • Tsunetsugu and Motome HE] suggested 
alternating d xy -d xz and d xy -d yz orbital chain structure, which is compatible with the 
space group IAx/a. On the other hand, Tchernyshyov [19] argued against the above 
model and proposed the space group IA\/amd which is associated with a complex 
linear combination of d xz and d yz orbitals. Regardless of the detailed picture, the one- 
dimensional antiferromagnetic chain is the essential ingredient to both models owing to 
the direct overlap of d xy orbitals [I]. This means that the orbital ordering imposed by 
the exchange couplings of V 3+ — Mn 2+ ions is different from the orbital pattern favored 
by the structural phase transition. This causes the instability of the antiferro-orbital 
configuration. Although frustration is relieved through the structural phase transition, 
it also paves the way for the competing orbital ordering of V 3+ ions to occur. 

Lastly, we present the external field (H) dependence of the NMR spectrum in figure 
6. We observe that each identified line at zero field does not shift at all with H. Instead, 
amplitude of each peak decreases or increases without shift in such a way that the first 
moment (Mi) or the center of gravity of the spectrum changes linearly with increasing 
H as shown in figure 6 (c). Moreover, new NMR lines appear as H increases. The slope 
of Mi versus H is initially 5.4 MHz/T but get bigger smoothly up to 3 T being fixed 
to 7 MHz/T. The value is still smaller than the gyromagnetic ratio 77V of V (11.193 
MHz/T). We conjecture that this behavior may be possible only with assumption of the 
discrete spin directions which is distributed spatially. The change of the slope below 3 
T may be due to the effect of the rotation of the domain or the similar stabilized process 
against the external field. 

Also we observed that the total intensity of the spectrum increases with increasing 
H up to 4 T above which it is saturated [see figure 6 (b)]. Typically, NMR signal 
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Figure 6. (Online color) (a) External field dependence of spectrum up to 5 T. (b) 
Total integrated intensity increases with external field saturating at ~ 4 T. (c) First 
moment versus external field. The initial slope of 5.4 MHz/T changes to 7 MHz/T at 
higher fields. The value must be compared to the expected value of 11.193 MHz/T. 

enhancement is associated with the domain walls which are swept out with H resulting 
in the disappearance of the enhancement effect. In this regard, the usual enhancement 
of the signal with large H suggests the persistence of the domain motions at least up to 
4 T. Since the domains of the long-range ordered phase are expected to be frozen-in at 
low fields, the enhanced high-field signal might be related to the domain motions of the 
short-range ordered state. Above we have ascribed the short-range order to the orbital 
ordering induced by the structural phase transition. Specifically, the V ions form the 
antiferromagnetic chain with the exchange coupling constant of an order of 1000 K [3]. 
In this situation, the short-range ordered component of the V spins is frustrated. Thus, 
in the studied field range the short-range ordered state undergoes a minor change in the 
spin directions accompanying the domain motions. To draw the more concrete picture 
on the puzzling behavior of the field dependence, a further study of the single crystal is 
needed. 

Within our study, we cannot determine the exact nature of the short-range ordered 
ground state. More direct determination has to come from a neutron diffraction study. 
However, it cannot be identified as a classical spin glass phase because its NMR signal 
intensity is comparable to that of the long-range ordered phase. It is possible that the 
RSG-like behavior is associated with the FEM domains and is a consequence of the 
freezing of fluctuations of the short-range ordered state. 
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4. Conclusion 

In conclusion, we examined the ground state of MnV204 through magnetization and 
NMR measurements. We observed the reentrant-spin-glass-like behavior on top of the 
long-range FEM order. Our study provides a clear-cut example of the inhomogeneous 
ground state in undoped compounds where spin interactions and orbital ordering related 
to a structural phase transition are two competing parameters. 
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